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I. INTRODUCTION 



The electromagnetic processes influenced by the Abelian axial anomaly [l|, |2| are of 
considerable theoretical interest. Among them are the transitions of the type 'y*{q) — >■ 
P^{pi)P^{p2)P~{p3), where 7* denotes a, generally, virtual (g^ 7^ 0) photon 7, stands 
for a charged and for a neutral meson from the pseudoscalar nonet, up to the strangeness 
conservation (so that P^ = 7i^,K^ and = tt^,?],?]'). These processes are supposedly in- 
fluenced by the, colloquially called, "box" axial anomaly, since on the microscopic level, 
the three pseudoscalar (P) mesons would couple to the photon through a four-vertex quark 
loop, like in Fig. [1] 

In the chiral limit (where m^r = 0) and the soft-point limit (of vanishing 4-momenta of 
external particles, Pj = = q), which is a reasonably realistic approximation at low energies 
at least for the lightest pseudoscalars - the pions, the anomaly analysis predicts jsj-lGj that 
the theoretical amplitude is exactly 

where e is the proton charge, A''^ the number of quark colors, and the pion decay constant 
U = (92.42 ± 0.33) MeV, whereby = (9.72 ± 0.09) GeV-^ 



7r+(pi) 




vr (233) 

FIG. 1: One of the box diagrams for the process 7* — t- tt+tt^tt". There are six different 
contributing graphs, obtained from Fig. [1] by the permutations of the vertices of the three 
different pions. The position of the u and d quark flavors on the internal lines, as well as 
Qu or quark charges in the quark-photon vertex, varies from graph to graph, depending 
on the position of the quark-pion vertices. The physical pion fields are tt^ = (tti =p z7r2)/-\/2 
and 71^ = TT^. Thus, in Eq. ([3]) one has vr = \/2{n^t^ + n^tJ) + ix^t^ where t± = ti ± it2 

(see text). 
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On the other hand, the experimental knowledge of the processes that should be influenced 
by the "box anomaly" is not at all satisfactory, being quite scant. For the 7* — > 7r+7r°7r~ 
processes, which should be best approximated by the anomaly prediction ([1]) since it in- 
volves only the lightest pseudoscalars, there is only one published experimental value for the 
amplitude at finite momentapj, i.e., the form factor F!^'^{pi,p2,P3). It was extracted from 
the cross-section measured |7| at Serpukhov in the transition 7r^7* — t- tt'^tt^ through the 
Primakoff effect, so that its value F^'^[expt) = 12.9 ± 0.9 ± 0.5 GeV~^ really corresponds to 
the average value of the form factor over the momentum range covered by the experiment. 
The 7r~ scattering on electrons at CERN SPS yielded the total cross section |8|consistent 
with the Serpukhov value. (It is maybe cautious to recall that the both results [71, [8| are in 
the strong disagreement with a rarely quoted analysis jsl of an old measurement |10| of this 
elusive process.) In the meantime, one still awaits the analysis of the measurements of this 



form factor performed at CEBAF 11 | 



Now, however, there are new hopes of more and better experimental knowledge of such 
processes, as new high-statistic data on the form factor for the 7r~7* — )■ 7r~7r° transition 



are expected soon from the COMPASS Primakoff experiments at CERN [12|, |13|. (Not 
only pion, but also kaon beams can be used in these experiments, so that also the reaction 
K-Y -> K-n° can be studied by the COMPASS collaboration.) 

Thus, experiments may finally confirm the relation ([2]) between the "box anomaly" pro- 
cesses and much better understood and measured "triangle anomaly" processes, notably 
the vr°(p) —J- 7(^:1)7(^2) decay into two real photons, kf = = k^. Namely, the pertinent 
chiral-limit and soft-point amplitudes = limm^^o ^7r(/^i, ^^2) and A^^ are related jsl-H] as 

Al' = eflAt\ (2) 



where, in the notation^ of Ref. IJ], T,r(^i, ^2) is the (unnormalized) vr'' — )■ 7*7* form factor. 

The axial anomaly, which dictates these results, occurs on the level of the quark substruc- 
ture of hadrons for the quark loops to which pseudoscalar mesons are coupled through an 
odd number of axial (A) vertices, while photons are coupled through vector (V) vertices. On 
the level of effective meson theories where quarks are completely integrated out, the effects 
of the axial anomaly are encoded in the Wess-Zumino-Witten (WZW) Lagrangian term 



15|, ll6|. Nevertheless, thanks to the Veltman-Sutherland theorem^, the anomalous ampli- 
tude A^ is also obtained successfully through the non-anomalous PVV triangle diagram, 
where the pion is coupled to the fermion loop through the pseudoscalar (P) vertex. In the 
simplest variant, this is achieved through the, basically, "Steinberger-type" [l9[ calculation 
supplemented by the quark-level Goldberger-Treiman (GT) relation gj^gg/Mq = l/Z^r con- 
necting the (constituent) quark mass parameter Mq and the vr-quark P-coupling strength 
g^qq with 

Such simple "free" constituent quark loop (CQL) calculations are surprisingly successful. 
While in the present context the most important is their exact reproduction of the "triangle" 



^ Except that here, because of F!^^ {pi,p2,P3) and related function depending on three independent mo- 
menta, it is for brevity not written explicitly, but understood implicitly, that the scalar functions depend 
only on the scalar combinations of momenta; e.g., r^(fci, fc2) = T7r(fcJ, for on-shell pions. 

^ Namely, this theorem dictates that in the chiral limit, the PVV tt" — >■ 77 amplitude is given exactly by 
the coefEcient of the anomaly term (see, e.g., 0, Q)- 
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and "box" anomalous amplitudes in ([T]) and ([2]), let us also recall that just the PW quark 
triangle amplitude leads to over 15 decay amplitudes in agreement with data to within 3% 



and not involving free parameters 18|, |20|, |2l| . Since "free" quarks here mean that there are 



no interactions between the effective constituent quarks in the loop, while they do couple to 
external fields, presently the photons and the pions tTq, the simplest CPT, Lorentz and 
SU(2) invariant effective Lagrangian encoding this is^ 

^cff = -q{^ -ieQ^ + + 2z^^g-,757r + ...)q, (3) 

where tt = ■nata, ta = Ta/2 and Pauli matrices, whereas Q = diag(|, — |) is the charge 
matrix of the quark iso-doublet q = (m, d)'^. The extension to SU(3) is obvious. The resulting 
CQL model calculation would be the same as, e.g., the lowest (one- loop) order calculation |2| 
in the quark-level linear a-model j22|,[23|. Hence the ellipsis in the Lagrangian (I3]) - to remind 
us that Eq. ([3]) also represents the lowest order terms pertinent for calculating photon-pion 
processes, from the a- model Lagrangian and from all chiral quark model Lagrangians (e.g., 
see [23]) containing the mass term with the quark-meson coupling of the form 

-M,q{UPL + U^PR)q, (4) 

where Pl,r = (1 ±75)/2. Namely, expanding 

[/(t) ^ exp[(-)z7r/2A] (5) 

to the lowest order in pion fields and invoking the GT relation, returns ([3]). 

In contrast to this simple CQL model, a more sophisticated approach to quark-hadron 



physics is provided by the Dyson-Schwinger (DS) approach (25l-l27l|. which has clear connec- 
tions with the underlying theory - QCD. Namely, this approach clearly shows how the light 
pseudoscalar mesons simultaneously appear both as quark-antiquark {qq) bound states and 
as Goldstone bosons of the dynamical chiral symmetry breaking (D^SB) of nonperturbative 
QCD, a unique feature among the bound-state approaches to mesons. Through D^SB in 
DS equation for quark propagators, dressed, momentum-dependent quark masses Mq(p^) 
are generated. They are in agreement with perturbative QCD for high momenta. However, 
thanks to D^SB, at low momenta they are of similar order of magnitude (and even tending 
to be somewhat higher) as typical constituent model mass parameters Mg. That is, Mq ~ ^ 
of the nucleon mass ~ ^ of the p-meson mass nip, or higher if the mass defect due to the 
binding of quarks is taken into account. This is true even in the chiral limit, i.e., for van- 
ishing masses of fundamental quarks, which underscores the nonperturbative character of 
D^SB. The DS approach thus provides a partial justification of this simple CQL model, and 
adds to the understanding of its aforementioned phenomenological success IJ]. Namely, 
although the CQL model obviously suffers even from a lack of some qualitatively essential 
features, notably confinement, the assumption is that below spurious qq thresholds, a more 
important role is played by D^SB (which generates large constituent quark masses, i.e., 
Mg ~ 300-500 MeV for u and d, the higher estimate being suggested by the DS approach; 



e.g., see Refs. |25l427l| and references therein). 



The metric is given by r/n = 1, 7722 = 1, ?733 = 1, '?oo = -1 
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The DS approach uses the solutions of Bethe-Salpeter equations for the pseudoscalar 
meson bound-state vertices instead of the point pseudoscalar couplings g-nqql^Ta of the CQL 
model ([3]). Then, both the anomalous amplitude and the connection ([2]) with the box 
anomaly amplitude are again (in the chiral and soft limit) reproduced exactly and 
analytically |28|, l29| and independently of details of dynamics, which is again unique among 
the bound-state approaches. 

The extension of these amplitudes from the chiral and soft-point limits to general form 
factor kinematics have often been studied; e.g., Ref. |30| used CQL (|3]) to study in this 
way the pre sently pertinent "box" amplitude. Present paper aims at continuing the study 
of Ref. 



30 



by examining the possibility of including also the vector mesons, primarily 
in the description of the "box" -anomalous transitions in a mixed quark- meson theory. As 
will become apparent below, this is a nontrivial and interesting theoretical issue in its own 
right, but there is also obvious phenomenological relevance in this context. For example, 
in the decays ?7, rj' — 27r7, the vector mesons turn out to be essential for reproducing 



experimental results in very different approaches such as 31| and 20[. Also, since in the 
process 7* — > Tr'^vr'^Tr— (further, 7* — Stt for short) one can depart strongly from low 
momenta, one may expect that the vector mesons will be important also here (e.g., see Ref. 
[3^). On the other hand, if a treatment of a process is phenomenologically successful thanks 
to the inclusion of vector mesons (as in the cases 77,77' — 27r7), an important question is 
whether anomalous processes are still described correctly in the low-energy (i.e., chiral and 
soft-point) limit. This problem has a somewhat lengthy history, of which only the presently 
necessary part will be reviewed in the beginning of the next section, which combines CQL 
with vector mesons and finds the resulting — )■ 27 and 7* Stt amplitudes, revealing 
superfluous contributions of the resonant graphs to 7* — Stt in this limit. In Sec. IIIIl we 
present a resolution of this problem. In Sec. IIVI we complete the calculation and discuss 
the results. We summarize in Sec. |Vl 



II. INCLUDING VECTOR MESONS 



A. Short history of vector mesons in 7* — t- Svr 



The Vector Meson Dominance (VMD) is certainly a reasonable approach to try because 
of its numerous empirical successes regarding electromagnetic interactions of hadrons (e.g., 



see 



33| for a review and references), although its basis in QCD has not been fully clarified 



yet. In a purely mesonic theory, Rudaz S^l assumed VMD where interaction with photons 



takes place only through p° and u mesons. For example, 27 would occur through 



7r° — p^uj and p° -> 7. u; — 7 [35 



Assuming the appropriate relationships between 
the pertinent coupling constants, he successfully reproduced the anomalous amplitude for 
tt" — )■ 27. Nevertheless, with the standard Kawarabayashi-Suzuki-Fayyazuddin-Riazuddin 
(KSFR) relation [36|, l37| between p-meson mass rup and fmn coupling (7p7r,r, 



9, 



pTTTT 



vr 



(6) 
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this VMD approach ( "pure" VMD in the following) would then give the amplitude A^^ too 
larg e by the factor of |, violating the axial anomaly relation ([2]). To avoid this, he advocated 
in agreement with Zinn- Justin and collaborators jssj) the KSFR relation revised by 



the factor ^ 34 



f within few 



However, the experimental values of gpnn, ''^p and f-^ strongly support 



Vo) the original one ([6]), which Rudaz finally adopted, also introducing 39[ the 
contact term for the direct u ^ Sir transition. Namely, for a favorable choice of its coupling 
strength, this term could contribute — | of the correct amplitude A^"^, finally enabling the 
VMD approach to reproduce 39|] the axial anomaly predictions ([1]) and ([2]). We denote this 
by "modified" VMD. Subsequently, Cohen showed 40| that the pertinent Ward-Takahashi 
identities (WI), first derived by Aviv and Zee joj, support the existence of such an extra 
contact term. 



B. CQL VMD models 



We want to examine whether anomalous processes like 7r° — 27 and 7* — )■ Sir can be 
properly described by a mixed model of constituent quarks and mesons, which, unlike the 
Lagrangian (|3]), would include not only pseudoscalar, but also vector mesons. Thus, the 
quark-meson-interaction part of the Lagrangian ([3]) gets enlarged to 



int 



(7) 



where'^ = p^ta + ojHq, = ro/2, tq = diagfl, 1). This quark-meson interaction is, for 
instance, used in the quark-loop approach of Refs. [18|, |20|, |21[ - a prominent example, since, 
as pointed out earlier, it describes many processes without involving free parameters. 




7r+(pi) 



FIG. 2: The form factor GpjrTr seen as a quark triangle diagram in the CQL model. 

With this theory at hand, one could compute the Gptttt form factor in the quark-loop 
model, where the p — t- tttt amplitude is represented by a quark triangle (Fig. [2]) , and its 
crossed mate, with VPP couplings. In terms of an effective Lagrangian, the pvrvr interaction 
would then be given by 

^p^^ = 2zG'p^^Tr([7r, {dp'K)]p^), (8) 

which implies that besides the "box" VPPP graphs in Fig. [H the process 7* — )■ Sir also 
receives contributions from the p-resonant triangle graphs such as the one in Fig. |3l 



We assumed the ideal uj-cf) mixing, as well as an U(2) symmetry for the interactions of quarks and vector 
mesons. 
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For GpTTTT we use the chiral and soft-point limit result by Hakioglu and Scadron [23 
limm-^^o^pmriQ = 0) = QpnTT = Qpqq = coust., cousisteut with the hypothesis of the VMD 
universality, which assumes the equality of all p couplings, namely the couplings to fermions 
(presently - quarks g), to pions tt, and to photons 7, i.e., 

gpqq = gpnn = 9 p^ = COUSt. = . (9) 

Electromagnetism is already present in the starting Lagrangian ([3]) through the direct 
quark-photon coupling ieqQjfLq, but to incorporate VMD in our theory, the coupling of 
the photon to neutral vector mesons must be added. One way to describe the interactions 
between photons and hadronic matter in the spirit of VMD is 

^VMDi = teqQj'-A^q - ^-F^^ (—p% + — ujp^ . (10) 

This version of VMD is often called VMDl. Here Fp^ = dp Ay — d^Ap, Up^ = dpOJy — d^Up, 
and also = dpp^ — 9^p° [s^. Ideal ui-(j) mixing gives (7^^ = 3gp^. VMDl can be 
transformed into the most popular representation of VMD where there is no direct quark- 
photon coupling. In the limit of univers ality (|9]), this standard ( "Sakurai's" ) representation 



of VMD, denoted VMD2, is equivalent 33|] to VMDl. The presently pertinent part of the 
VMD2 Lagrangian reads^ 

^vmd2 = ^AMpJ + -u;J. (11) 

Our starting CQL Lagrangian ([3]) is thus finally augmented to include vector mesons (further 
denoted as the CQL- VMD model): 

^cff = ^int + ^VMD2 + ^p.n " + + ...)q . (12) 




FIG. 3: One of the "resonant" triangle diagrams for the process 
7(g) — )■ TT~^ {pi)'n'^ {p2)'n'~ {ps) , where the two pions with momenta pi and p2 are obtained by 
the decay of the intermediate p meson transferring the momentum Q = Pi + P2- Two more 
analogous graphs are obtained by replacing p+ — >■ tt+tt" and vr" by, respectively, p° — )■ tt+tt" 
and 7r°, or by p~ 7r~7r° and tt"*". Also, each of the three graphs has its crossed graph. 



^ Where we omit the kinetic term and the fictitious "photon mass term" which arises (e.g., see 
the gauge invariance of tlie theory, which is manifest in VMDl ([TU]). 



331) from 
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The ellipsis in Eq. (fT2|) again serve to remind one that the terms not pertinent for photon- 
pion processes are not shown. (E.g., if the starting Lagrangian was the sigma-model 
one, we would have terms containing scalar, a-mesons also in Eq. flT2|) . but they do not 
contribute to the presently interesting 7* — > Sir and 27. That is, the contribution of 

the graphs analogous to the resonant triangles, but with intermediate a — tt+tt" instead of 
p° — TT+TT", vanishes due to parity conservation.) These terms are, in general, different for 
different theories, but since they are presently irrelevant, our conclusions will be the same 
for many different models, from the a- model to the chiral quark models 2J], after the p and 
u mesons are introduced to them. In that sense, CQL-VMD denotes not just one, but the 
whole class of models. 

In contrast to ^ymdi there is no direct quark-photon coupling in the VMD2 picture 
( IT2I) . Thus, in VMD2, the six box graphs exemplified by Fig. [T] do not contribute to 7* — Svr 
in the precise form depicted in Fig. [1] Instead, they are modified so that the photon first 
couples to the intermediate cj-meson which in turn couples to quarks. That is, the quark- 
photon coupling in Fig. [1] is replaced by the photon coupling gi^^ to u propagating to its 
qq vertex. There is no need to spend space on re-drawing Fig. [1] to depict this insertion 
of uj since it modifies the photon coupling in the same way in all graphs, and is illustrated 
in Fig. |3l (Note that in VMDl, there are both the graphs with the direct quark-photon 
coupling and their partners with the cu-insertion, but since the momentum dependence is 
different than in VMD2, the sum of these graphs in VMDl yields the same results as VMD2.) 

The box and resonant triangle graphs in which the photon is coupled to p° (instead of to 
u), all vanish, as must be due to G-parity conservation. 



C. 7r° 27 through CQL VMD models 

Since we choose to work with VMD2, the outgoing photons in this process are created 
only through the mediation with the and u. In our model, these, in turn, come from 
a "triangle" PVV quark loop. If the photons are on-shell there is a complete cancellation 
of the u and p propagators with in the VMD coupling. In the same manner, and by 
using universality, the gpqg and g^^gq couplings get canceled with the gp^ and g^j-y, respectively. 
Then, our Lagrangian f lT2|) leads to the same 7r°(p) — )■ 7(A;i)7(A;2) amplitude as the standard 
quark-triangle-loop calculation (e.g., see jssj), namely 



= h^e:{ki,ai)k2pelik2,a2)e^''''^T^{h,k2), (13) 
where fcj and o"j are the momentum and polarization of '~f{ki), and 

TAkuk2) = ^^^o(fci>^2). (14) 



Here Co = (2!M^/z7r^)Co, where Cq is the standard 't Hooft-Veltman 4l| scalar three-point 
function. The limit <^ Mg, together with the GT relation, reproduces the analytical 
result 

lim = ^ _ = Al^ (15) 
" 127r2 " ^ ^ 

for the anomalous chiral decay into two real photons, kf = 0, /c| = 0. 
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D. The problem with 7* Svr in CQL VMD models 

Previous short calculation served as a consistency check aimed at reproducing the cor- 
rect low energy limit even when vector mesons are introduced. Here we perform a similar 
calculation for 7*(g) '^'^ {pi)t^^ {P2)'^~ (ps) , and obtain the amplitude 



= e^(q,a)pi,p2,p3A£^"'"i^f (pi,P2,P3), (16) 
In this notation, the contribution from the resonant triangles reads 



^2 ^2 

FnPuP2,P3) 



1 eiVc 97rqq 9 pqq 

2 Gtt^ Mq ml ml + 



X 



2 2 2 

771 ^ mi ~ mi ~ 

-Y^Co{pi,P2) + ^^Co{p2,P3) + ^^Co{puPs) . (17) 

Lm^ — s m^ — t m^ — u 



where the Mandelstam variables are defined^ by s = —{P1+P2Y, t = ^(^2+^3)^, u 
— {pi + PsY- The form factor from the box graphs of Fig. [1], is 



FoiPuP2,P3) = 1^3nqq9lqq^;r}^ 

' ' 3 Qtt"^ Mq m^p m^p + 



X 



Do{pi,P2,P3) + Do{pi,p3,p2) + Do{p2,Pi,P3) (18) 



where Dq = (3!M^/i7r^)i5o and Dq is the 't Hooft-Veltman scalar four-point function 41 . 
The total amplitude would be F^^ = F'^^ + F^. 

In the soft-point and chiral limit (i.e. pi,P2,P3 — > 0), where Co, -Do — > 1, and with the 
usage of GT and KSFR relations, 

pres 3 eAT, 1 _ 3 3^ F ^ eAT, 1 

-"2 12^;^- 2"^^' ^<^^T2^:^ - ' ^^^^ 

which would mean a total of ^^^^^(0,0,0) ^ | f^jj = | A^'^, which is by the factor 5/2 
bigger than the correct anomalous amplitude. Notice that if we had included only the 
resonant triangles we would have a result that is off by 3/2; these are the very same 3/2 
that we mentioned after Eq. ([6]) as the reason for Rudaz adding the u — Sn contact term 
into the previously "pure" VMD description. It appears that the above Lagrangian ([7]), 
( ITTl) leads to an inconsistency; if we want to calculate anomalous processes with more than 
one pseudoscalar, it is not legitimate to just add vector mesons a la Sakurai to a mixed 
meson-CQL. It is easy to check by explicit calculation that the same problem persists (as 
expected) if VMDl approach is used instead. 



Note that this definition is different from the eonvention in the Serpukhov pap er [7| and, e.g., Rcf. |42|. 



Our notation is closest to that of the proposal of the CEB AF experiment [ll| , where the outgoing pion 
pair is tt'^tt'^. Thus the choice s = — (pi +^2)^ = ^{Pir+ + Pti")'^^ while the squared invariant mass of the 
pion pair outgoing in the Serpukhov experiment, 7r~7r°, is t = —{p2 +^3)^ = --{Ptt° + P-k-Y (see Fig. [T]). 



9 



III. A RESOLUTION THROUGH WEINBERG TOMOZAWA INTERACTION 



In the present context of strong interactions and hadrons, which, after all, have substruc- 
ture, one may think of introducing a form factor [Q"^) for the transition from the vector 
quark vertex to 2ti instead of the corresponding part of the resonant triangle graphs (Fig. 
[3]), which makes the latter troublesome. Such a form factor would be constrained by the 
known anomalous behavior of the total form factor F^'^(pi,p2,P3) in the soft point. How- 
ever, introducing it by hand and ad hoc, without an insight into the underlying dynamics, 
would not be a satisfactory way of removing the superfluous contributions. To understand 
which modifications of the CQL-VMD approach to make, we seek guidance from a more 
fundamental, substructure level. 




FIG. 4: One of the graphs for the 7* — )■ Svr non-resonant triangle contribution depicting a 
7r+7r° pair coming out from the qqirir Weinberg- Tomozawa (WT) interaction that is needed 
for the appropriate subtraction. The two analogous graphs are obtained by exchanging the 
pion pair by tt+tt" and tt^tt". There are two of each of these graphs. 



A. Insight from the Dyson- Schwinger approach 

In the Introduction we already mentioned the DS approach to QCD. Although there is no 
full derivation of CQL or VMD from QCD as the underlying theory, much of their features 
can be reproduced and understood in the process of describing pseudoscalar and vector 
mesons and their interactions in the DS approach, which is free from the above problem of 
the superfluous contributions | f lT9|) of the resonant graphs to the anomalous amplitude. 

The DS approach employs dynamically dressed quark propagators S{k) = [i^A^k"^) + 
obtained by solving the "gap" DS equation, so that the momentum-dependent 
mass function MqihP') = B {k'^) / A^k"^) takes place of the simple constant constituent mass 
Mg. Also, the solutions of the Bethe-Salpeter (BS) equations for the pseudoscalar meson 
bound-state vertices replace the point pseudoscalar couplings g-^qqlb of the CQL model ([3]), 
and a vector Wl-preserving dressed quark-photon vertex T^{k,k') is in place of the bare 
quark-photon vertex 7,,,. (In the vast majority of phenomenological applications this is the 



Ball-Chiu Ansatz 25|, |26|.) This procedure defines the generalized impulse approximation 



(GIA). As we already mentioned, this reproduces the A'^"' and A^^ low energy theorems. 



10 



On the other hand, the momentum dependence, i.e., the growth of the form factor 
F!^'"{Pi,P2,P3) from its soft-point hmit Ai^'^ turned out to be so slow (after all six per- 
mutations of the graph in Fig. [T]were taken into account properly j29|), that the Serpukhov 
data point could not possibly be explained in that approach. (Also the related processes 
?7, rj' — > 27r7 could hardly be reconciled with such a weak momentum dependence 29|| of the 
7* —J- Svr form factor.) This indicated that the DS approach should include the contribu- 
tions from vector mesons - or rather, in the DS context, their microscopic, gg-substructure 



equivalent. To this end, Ref. |43j went beyond GIA in its treatment of the box graph, 
inserting and summing the infinite set of gluon ladder exchanges in the s, t and w-channel 
of the box graph. The resulting momentum dependence of F^'^{pi,p2,P3) agrees well with 
the VMD behavior; the inclusion of these gluon diagrams beyond GIA successfully produced 
p-meson-like intermediate states in the two-pion channels. Thus, their contributions cor- 
respond to our resonant triangle graphs, but with the crucial difference that in the DS case 
these contributions vanish as one approaches the soft limit of vanishing momenta, so that 
the correct anomalous amplitude is obtained, unlike Eq. (fT9i) . 

This favorable behavior can be understood on the basis of Maris and Tandy's DS results on 



dressed vector vertices [4J] . In the DS approach one does not have elementary meson fields, 
but one explicitly constructs physical, on-mass-shell mesons as qq bound states which are 
eigenstates of mass. Thus, the mesons are not well-defined away from their mass poles, nor 
are their couplings (such as gpnn, Qp-y and gpqq). Therefore, at the level of meson substructure, 
the issue of the p-meson-like intermediate states, including the resonant p contributions to 
the two-pion channels, should be addressed within the dressed quark vector vertex r^(A;, k'), 
which couples not only to photons but also to the pion vector current Tr([7r, (9^7r)]ta ). 

Using the same DS dynamical model^ for the quark-gluon interactions as later Ref. jisj , 
Maris and Tandy 4J] solved this inhomogeneous BS equation for the dressed quark vector 
vertex r^(fc, /c'). In this way they essentially reconciled the VMD picture with the QCD 
picture of a photon coupled to quarks in a qq bound state; for the present paper their most 
important result was that the inhomogeneous BS equation generated both resonance and 
nonresonance contributions to the full (model) vector vertex Tp{k,k'), which, unlike the 
BC vertex, contains timelike vector-meson pole (at the model value = —my = —m? = 



-m: 



2 



-0.55 GeV^) in the part of the vertex transverse to Qp. While this part (the resonant 
part) of the vector vertex BS solution is significantly enhanced over the BC Ansatz, it also 
vanishes as — )• 0. This is explicitly shown in Eqs. (30) and (34) of Ref. ji^ ]. 



B. Subtraction of the Weinberg- Tomozawa interaction 

In terms of meson degrees of freedom, this means that the resonant contribution from the 
intermediate vector meson with = is absent. This is reasonable as it would correspond 
to a constant meson propagator l/my, and this in turn corresponds to a point interaction 
in the coordinate space. In the present case, it would correspond to a p-meson propagating 
zero distance from its qq vertex before turning into two pions, i.e., to a qq vertex producing 



^ This model is the most widely used one in the phcnomenological branch of DS studies - see DS approach 
reviews such as Refs. 25l-l27| 
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two pions immediately, since here is either s, t, or u, depending on the two-pion channel 
coupled to p. Note that this unphysical situation is quite different from the situation when 
such "intermediate but non-propagating" p turns into a 7: this just means that also VMD2 
incorporates implicitly the situation when the photon couples to quarks immediately and 
directly (for example, as in the case of 7r° — 27, where the factors of rup and couplings, 
except e, cancel). Contrary to that, two pions can come from a vector qq vertex only via a 
truly propagating intermediate p meson (with transferred momentum Q"^ ^ 0), while there 
is no direct conversion of quarks and antiquarks into the two-pion vector current. 

This is the reason why our previous calculation has led to the spurious soft-point contri- 
bution I A^'^ (HM : in the resonant graphs, the part 

Ql'^taQ intermediate — )■ eabcT^bd^T^c 
yields simply^ the "vector qq — )■ 2n form factor" 

1 Tfl^ 1 

(Q ) ^ 9pqq 2 ^P'^'^ ^ o72 7l2~i 2 ' 

^ + m^p ' 2/^ + 

which is, nevertheless, wrong as it stands because it contains the contribution of the in- 
termediate p with = 0. The correct vector qq — )■ 27r form factor can be obtained by 
subtracting this contribution: 

^^(g^) = ^;(g^) - ^;(0) cx gp,, gp^^ - gp,, gp^^ = ^ ^^^^ , (21) 

i.e., the resonant contributions depend on essentially as in the DS substructure consid- 



erations such as j43| and j4J| (see esp. Eqs. (30) and (34)). 

One may visualize the removal of the point qq — > 27r interaction (non-propagating, = 
p) leading to Eq. ( l2T]) as the subtraction of the point-interaction triangle graphs (such as 
Fig. H]) from the corresponding resonant triangles (e.g.. Fig. [3]). In terms of formulas, this 
subtraction corresponds to including the following vrvrgg point coupling into the effective 
pTTTT part like this: 

2igp^^Ti{[n, (9^7r)]p'^) 22Tr([7r, [dpn)]{gp^^p^' - ig^^q^J^)) (22) 

where = taJj^, J^^ = qi^taq, and g-^^qq is fixed precisely in a way to respect the 7* — )■ 37r 
low energy theorem, i.e., in a way that it cancels the resonance part in the soft-point limit 
completely, yielding g-^-^qq = Written in this form, it turns out to be nothing else 

TV 



but the quark-level Weinberg- Tomozawa (WT) interaction [45|, |46|. However, it must be 
understood that Eq. ( l22l) indicates only the subtraction of the resonant graphs and does 
not mean adding a new nnqq interaction term to the Lagrangian ( fT2l) . If one would try 



this, the KSFR relation would be spoiled in the same way as when the analogous two-pion 



interaction with nucleons is added to the VMD-nucleon Lagrangian 40 . 



- thanks to vanishing of the transverse part of p propagators when contracted with Levi-Civitas from 
the traces of PVV triangles. The second equahty is from the KSFR relation. 
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This subtraction, that needs to be included only when dealing with P W triangles coupled 
to vector mesons decaying to two pseudoscalars, completes the definition of our constituent 
quark model coupled to pseudoscalar and vector mesons. With it, the combination of the 
resonant triangles (e.g.. Fig. [3]) and the "subtraction triangles" with the added point vrvrgg 
interaction (e.g.. Fig. S]) yields the behavior in accord with the Abelian axial anomaly 
of QCD, as shown in detail in the beginning of the next section, where we complete the 
calculation of the 7* — ?■ Stt form factor. 



IV. RESULTS AND DISCUSSION 



A. Completing the calculation 



When we include the non-resonant, WT triangles (Fig. H]) in the amplitude, the complete 
triangle form-factor Fa = + F^"^ becomes 



i^A(pi,P2,P3) 



1 eNcgnqgd, 



qq ilpqq 



2 Gvr^ Mq m? m? + 



s 



— s 



Co{pi,P2) + 



— t 



u 



nip — u 



CoipuP3) . (23) 



With the usage of GT and KSFR relations, the constant prefactor becomes 



eNcQ^qq gl _ eNc 



67r2 Mq ml Uw^f^ 



A: 



The total amplitude is F^'^ = F^ + Fq. In the soft-point limit (i.e. pi,p2,pz 0), where 
box contribution 



F^^^ + -F^^ —J- and the total is in this limit given by the pseudoscalar 



^7^(Pl,P2,P3 ^ 0) ^ Fo(pi,P2,P3 ^ 0) ^ 



eN, 1 
12vr2 /3 



A, 



Stt 



It is interesting to note that for real photons (g^ = 0), and by "squeezing" the quark 
triangles and boxes to points, the total form factor F^'^ = Fo + Fa reads 



Svr 



+ 



t 



mi 



m^ — t 



+ 



u 



— u 



(24) 



where we used the GT and KSRF relation. This is precisely the Terent'ev phenomenological 
form 0, 13 for (5 = and Cp = 1/2 (in other words, a "modified" VMD result H Isoj 



The following Figs. [5H7|show the normalized form factor F^'^ = F'^'" j . For constituent 
u and d masses the typical estimates are around Mq ~ m^/S ^ 330 MeV and Mq ~ mp/2 f« 
385 MeV, but Figs. [5H7] also show results starting from Mq = 300 MeV and going up to the 
DS scale 28|] A = 565.69 MeV. The 't Hooft-Veltman integrals, Cq and Dq, were calculated 
numerically for t he case where the photon can be taken on-shell ~ (pertinent in all 



experiments (7|], [ll|], |l3|), and for 



13 




FIG. 5: (color online) For various values of the mass parameter Mg, our numerically 
calculated form factors F^'^(s,t, are shown as functions of t = — (p2 +^3)^ 
= —{PnO +Ptt-Y, i-e., of the invariant mass of the outgoing pion pair vr^vr" in the 
Serpukhov and COMPASS experiments. There, all pions are on shell, and for definiteness 
we fix u = m^. The curves belonging to the lower stripe (blue online) are only the quark 
box contribution (i.e., the lower stripe represent predictions of the simple CQL approach 

for various Mq). The upper stripe (with black curves) represents the corresponding 
predictions from the presently pertinent CQL-VMD approach; that is, the black curves 
represent predictions which include the p-resonant triangle loops with the subtracted 
Weinberg- Tomozawa interaction. In the both stripes, the dotted curves correspond to 
Mq = 300 MeV, the dashed ones to Mq = 330 MeV and the solid ones to Mq = 360 MeV. 
The exhibited data point j3] really corresponds to the average value of the form factor over 
the momentum range covered by the experiment (between the two-pion threshold and 10 
m^), which measured the total cross-section - see the discussion in Subsec. IIVBI 
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explored at CEBAF will be mostly in the s channel, and the amplitudes themselves 
have a weak dependence on the virtuality of vr", so we take ~ m^. Now s + t + u = 
ml, and we also fix t = — m^. 

Quarks are not confined in our model, so there are possible spurious contributions to the 




4.0 6.0 8.0 10.0 12.0 14.0 16.0 



t [ml ] 



FIG. 6: (color online) The form factors F!^^{s,t,u) from various approaches are depicted as 
functions of t, the invariant mass of the outgoing pion pair vr°(p2), T^'iPs), as in the 
Serpukhov and COMPASS experiments, where all pions are on shell. We fix m = for 
definiteness. The upper shaded stripe covers the results of our CQL-VMD approach for 

constituent masses between Mq = 360 MeV (corresponding to the solid black curve 
marking the upper edge of that stripe) and Mq equal to the DS scale A = 565.69 MeV of 
Ref. |29| (corresponding to the lower edge of that stripe). The lower shaded stripe (blue 
online) covers the results of the "pure" CQL model 30|] for the same Mq interval. That is, 
the (blue) dashed curve depicts the CQL model form factor for Mg = 360 MeV, while the 

lower edge of that stripe is the very slowly varying CQL form factor for the high 
Mq = A = 565.69 MeV, the DS scale of Ref. [29|. A comparison is made with results of the 



"modified" VMD [32, l39( (green dotted curve) and of DS (in GIA) [29( (red dash-dotted 
curve). Again, the exhibited Serpukhov point [3] is actually the average value extracted 
from the total cross-section - see Subsec. IIVBI for the comparison with experiment. 
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amplitude from the qq channel in the box as well as in the triangle if any of the Mandel- 
stam variables s, t, u, is bigger than 4M|. For Mq varied in the range of some 300 — 500 
MeV, the qq thresholds are at 600 — 1000 MeV. Even the lower value, 600 MeV, is beyond 
CEBAF upper bound [llj: y/s^ = Am^ = 554 MeV (for = 138.5 MeV). The vector 



mesons are especially important at COMPASS, where the proposed momentum range to be 



covered should go well above the vr vr^ threshold, up to the p-peak 47| . For such momenta 



constituent quark masses Mg employed in the CQL-VMD approach should be Mq > mp/2, 
say around 400 MeV, in order to avoid spurious qq thresholds. 

The general behavior of the amplitudes can be understood through three essential factors. 
The first is the presence of the p resonance, which is the dominant cause of the increase of 
the amplitude in the present approach. 

The second factor is the characteristic mass scale of a given model: in the simple CQL 
(e.g., Ref. |30|) and the present CQL-VMD models, this scale is simply the constituent 
mass parameter Mq, which is typically not much higher than ~ |M, - ^ 
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FIG. 7: (color online) Same as previous Fig. [6l but for the CEBAF kinematics, where the 

two outgoing and on-shell pions are vr^ and 7r°. The form factor is thus given as the 
function of their invariant mass squared, the s- variable. For the off-shell pion, vr", we use 

p\ = m^. We also fix t = —m'^. 
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models, their characteristic scales are also related to their dynamically generated momentum- 
dependent constituent masses a nd are relatively hi^ compared to typical Mg (e.g., A = 



565.69 MeV in the DS model 28| used in Ref. 29|). In the CQL model, momentum 
dependences are stronger for smaller values of Mg, while for Mg = A = 565.69 MeV it 
is even slightly weaker than in the DS approach with this scale A 29|. (See Figs. O- 
[3 and Ref. 13.) That larger characteristic mass scales suppress more the momentum 
dependence of -^^'^(^'1,^25^3)5 is manifest in its power series expansions in CQL [30| and 



DS 29| papers, where vector mesons were not included. Nevertheless, it is sufficient to 
note that the contributions of the quark loops are suppressed if the quark propagators are 
suppressed by large masses in their denominators. Thus, understandably, the present CQL- 
VMD approach also exhibits weaker momentum dependences for for larger M^'s, although 
now the p-resonance of course dramatically boosts this dependence overall. 

The third factor is the symmetry of the 7* — )■ Sir amplitude under the interchange of 
the external momenta p,. It was shown in earlier CQL 30|] and DS 29|] approaches without 
vector mesons, most clearly in the aforementioned expansions of F^'^jy^ , p^, y?) in powers 



of the momenta Pi (divided by an appropriate mass scale). Ref. [29| clarified how due 
to this symmetry, the contribution of the terms of the second order in momenta [(9(p^)], 
is in fact a small constant (of the order of m^) up to the virtuality of the third pion. 
Therefore, the main contribution, dominating the s, t, w-dependence for momenta smaller 
than some characteristic model mass scale, comes from 0{p^) and not 0{p'^). This gives the 
parabolic shape to the curves displaying form factors as functions of Mandelstam variables 



in Refs. 29|, |30|| and here. In conjunction with the large mass scale A, this also causes the 



weak momentum dependence of F^'" {pi,p2,pz) found [29|] in the DS approach using GIA, 
motivating Ref. 43|] to go beyond GIA for this process. 

We also show the results for F^'" obtained from "modified" VMD (which is basically our 
Eq. (!24l) . see also 39|, and 3^ and references within), and DS approach in GIA 29 . 



B. Comparison with experiment 

We should first note that the Sepukhov experiment 0] did not, in fact, measure the 
presently pertinent form factor F!^'"{pi,p2,P3), but the Primakoff total cross section atot- 
The latter is thus the experimental quantity which is the safest to compare with various 
theoretical predictions, which we do in Fig. [HI The measurements on various targets (with 
eZ being the nucleus charge) yielded atot/Z'^ = 1.63 ±0.23 ±0.13 nb [3], which is represented 
by the grey area in Fig. [8], where it is compared with the theoretical predictions of VMD 
and, for various constituent quark masses Mg between 300 MeV and 400 MeV, of the CQL 
and CQL- VMD approaches. The former is not compatible with the (admittedly scarce) 
experimental data, but its vector-meson extension, the CQL- VMD approach, is. 

This relationship between the CQL and CQL-VMD approaches would seem to indicate 
that agreement with experiment mandates the enhancement due to vector mesons and VMD 
already at the momentum scales where the Serpukhov data were gathered. However, such 



conclusion would be too rash, since Ametller et al. [43] showed that, in this kinematical 
domain, chiral perturbation theory (ChPTh) describes the measured 7r~7 — vt'^tt" process 
well, after the one-photon-exchange electromagnetic corrections are included (in their t- 
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channel, which is the w-channel in our conventions). This brought their theoretical prediction 
for the cross-section to a^^/Z'^ = 1.33 ± 0.03 nb, in agreement with experiment. They also 



pointed out the importance of the one-loop [49|] and two-loop 50| chiral-perturbation-theory 
contributions, especially for extracting from the experimental cross section the form factor 
value (call it Fq^^p^ = A^'^^expt)) which would correspond to the unphysical soft point 



t 



u 



0). Table I in Ref. [42 



reviews how including 1- and 2-chiral-loop and 
electromagnetic corrections gradually bring about F^^^^^ = 10.7 ±1.2 GeV"'^, consistent 
with theory, Eq. ([1]). 

These electromagnetic- and loop-corrected ChPTh results indicate that at Serpukhov 
energies VMD is not yet needed for agreement with experiment. Thus, the fact that in 
the COMPASS measurements of 7* Sir, the momenta to be covered should surpass the 
Serpukhov range and approach the p-peak 47| in the vicinity of which VMD dominates, 
remains the strongest motivation to combine the simple CQL approach with vector mesons 
and VMD. 




300.0 



320.0 



340.0 360.0 

M_ [MeV] 
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FIG. 8: (color online) Total cross sections for various Mg predicted by CQM is depicted by 
the lower solid curve, while CQM- VMD yielded the upper solid curve. The (constant) 
value predicted by VMD alone is depicted by the dotted line, and the ChPTh prediction 

is given by the 



(electromagnetically and chiral-loop-corrected) of 1.33 ± 0.03 nb [42 

dashed line. The experimental cross-section is denoted by the grey area 
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V. SUMMARY 



The 7* — Stt form factor, presently being measured with hi gh s tatistics at CERN by the 



COMPASS collaboration through the Primakoff experiments [12|, ll3|, has been computed 
in the present paper using the simple free constituent quark loop model extended by vector 
mesons. This extension turned out to present a problem for the transitions connecting 
one vector and three pseudoscalar particles, because the box graph, which saturates the 
anomalous amplitude A^^ (in the chiral and soft limit), is then supplemented by the p- 
resonant triangle graphs yielding the superfluous contribution of f^^'^- The same problem 
appears in a; — )■ Svr and r], r]' — )■ 2%^ decays, where the contribution of the p-resonant triangles 



successfully reproduces the empirical decay widths 18|, |20|, but the chiral and soft-point 
limit of the pertinent amplitudes are then in conflict with the low-energy theorems 

This problem is cured in the present CQL-VMD approach by removing the spurious 
contribution of the intermediate but non-propagating p-meson. Thanks to this, our model 
reproduces correctly the anomalous 7* — )■ Sir chiral- and soft-limit amplitude A^^ while 
including the effects of the vector mesons at higher momentum scales. 

The appropriate subtraction (|2T1) . of the form of the WT interaction 45|,|46|], was inspired 
by the insights obtained on a more fundamental, microscopic level through DS approach 



28, 29, 43, 44 



The question then arises why not use this, more fundamental QCD-based 
approach, to calculate the 7* — ?■ Sir form factor, instead of the present simplified approach of 
constituent quarks plus VMD. In fact, this was done a decade ago in GIA, but the momentum 



dependence of the resulting form factor F!^'^{s,t,u) is then very slow 29|], for the reasons 



explained in detail in the subsection IIV Al but also because VMD effects are lacking. Ref. 
431] thus endeavored to reproduce VMD effects working from a microscopic level, in DS 
approach, but found that this problem then requires going beyond GIA, making the task 
so intractable in spite of many model simplifications, that only the results for symmetric 
kinematics (with at least two Mandelstam variables equal) were given 43|. At this point, 
this DS approach (beyond GIA and reproducing VMD effects) seems hardly tractable for 
general kinematics, including those of COMPASS and CEBAF. Hence there is a need for 
related, more simplified models like the present CQL approach extended by vector mesons, 
called the CQL-VMD approach, where the required features are put in by hand under the 
guidance from phenomenology, WI and more microscopic approaches. Thus we may consider 
the CQL-VMD approach as mimicking the more microscopic DS approach 43[ beyond GIA: 
the relationship between the box graphs in the respective approaches is obvious, and the 
(resonant-subtracted) triangle contributions in the CQL-VMD approach may be regarded 
as mimicking the VMD effects reproduced in the DS approach beyond GIA by inserting and 
summing up the infinite set of gluon ladder exchanges in the s, t and u-channel of the box 



graph 43[. However, in contrast to the very demanding and difficult-to-use DS approach 
beyond GIA, the CQL-VMD approach can be easily used for any kinematics that are below 
the spurious quark thresholds. 

The DS approach beyond GIA, like in Ref. 43|] but for general kinematics, is a very 
difficult task which has to be relegated to a future work. The same holds for calculating 
electrornagnetic corrections, since in contrast to the corresponding calculation of AmetUer 



et al. 42], in the present framework the contributing diagrams contain the momentum- 



dependent form-factor F^^{s, t, u) in the 737r vertex. It makes this task much more difficult. 
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but it is obviously necessary in order to reach the next level of refinement in the present ap- 
proach. A much more straightforward future work in this direction will include a CQL-VMD 
calculation of the reaction K~'~f* — K~7i^, which can also be measured by the COMPASS 
collaboration 12|, |l3|. We will also test the present approach by applying it to numerous 



meson decays currently studied experimentally by WASA at COSY with high precision and 



statistics 



5ll . |52| , such as r], i]' . 
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